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ABSTRACT 

Direct numerical simulations of compressible, homogeneous, turbulent shear flows are 
used to evaluate Reynolds stress models. Three pressure-strain models, which are either 
linear, quadratic, or cubic in the anisotropy tensor, are considered. Dilatational dissipation 
and pressure-dilatation models are inserted into the Reynolds stress closure. Results show 
that variable-density extensions of incompressible pressure-strain correlation models do not 
correctly capture the compressibility effects seen in the direct simulations. In particular, the 
increase in the anisotropy of normal stresses and the reduction in the shear stress are not 
reproduced by any of the models. Also, the use of the incompressible form of the dissipation- 
rate equation to determine the solenoidal part of the dissipation is found to be questionable. 



1. INTRODUCTION 


A resurgence of interest in hypersonics has emerged, which is driven by advanced new 
applications such as high-speed civil transport aircraft, supersonic combustion ramjets and 
transatmospheric vehicles. These new applications bring into prominence some critical items 
that have limited the effectiveness of computational fluid dynamics codes used as tools for 
hypersonic system design. Chief among these is compressible turbulence modeling. Several 
turbulence models of varying degrees of complexity have been developed that range from the 
simplest algebraic or zero-equation model to the full Reynolds stress closure. Most of these 
models are simple extensions of their incompressible counterparts, where compressibility 
effects are incorporated into the models through changes in the mean density. However, 
many studies have shown that this type of model is unable to reproduce some features of 
flows that depend on compressibility, such as the reduction in the spreading rate of the 
compressible mixing layer as the Mach number increases. Thus, a better understanding of 
the effects of compressibility on flow turbulence is needed to improve current turbulence 
models. 

Recently, direct numerical simulation (DNS) of compressible, homogeneous turbulent 
shear flow by Blaisdell (1990) has shown that the growth of turbulent kinetic energy decreases 
as the Mach number increases. The reduction is due to two compressibility terms: the 
dilatational dissipation and pressure-dilatation correlation, which explicitly appear in the 
turbulent kinetic energy and mean temperature equations. Zeman (1990, 1991) and Sarkar 
(1991, 1992) have modeled the additional terms; the inclusion of the two terms in two- 
equation turbulence models leads to a significant improvement in predicting the reduction 
in spreading rate with increasing Mach number. 

In this paper, an assessment of Reynolds stress models in predicting compressible ho- 
mogeneous shear flow is conducted with the DNS of Blaisdell (1990). Three pressure-strain 
correlation models, which are variable-density extensions of their incompressible counter- 
parts, are considered. These models are either linear, quadratic, or cubic in the anisotropy 
stress tensor. The dilatational dissipation and pressure-dilatation correlation, which are the 
only explicit compressibility terms, are inserted into the Reynolds stress closure. Particular 
attention will be paid to the ability of each turbulence model to predict equilibrium states 
accurately. 

2. MATHEMATICAL FORMULATION 

In compressible, turbulent flows, two averaging techniques are commonly used to define 
the mean and fluctuating parts of a turbulent variable. Usually, conventional Reynolds 
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averages are used for the pressure p and the density />; Favre averages are used foi the 
velocity u, and the temperature T. Thus, any dependent variable / can be decomposed into 
mean and fluctuating parts in two ways: 

/ = / + /', / = / + /" (>) 


where the overbar represents the Reynolds average, the tilde denotes the Favre average, and 
the primes and double primes are, repectively, the deviations from the Reynolds average and 
the mass-weighted average. The Favre average / is a density-weighted Reynolds average 



( 2 ) 


We will consider the problem of compressible, homogeneous shear flow. In this problem, 
an initially decaying, compressible turbulence is subjected to a uniform shear S with the 
corresponding mean velocity gradients 


Uij = SS- il S 2 j 


(3) 


The flow is assumed to be an ideal gas, which satisfies the equation of state p - pRT. The 
evolution equation of the mean temperature, derived from the energy equation, is given by 


c v 


DT 

~Dt 


= £ - 


(4) 


where e is the total dissipation rate, iv d = p'd'/p is the specific pressure-dilatation, and c„ is 
the specific heat at constant volume. 

For a compressible homogeneous shear flow at high Reynolds numbers, the favre-a\eraged 
Reynolds stress tensor T tJ = u" u" is a solution of the transport equation 


^ = Pij + n.j - | eSij + (5) 

where n tJ is the deviatoric part of the pressure-strain correlation and Pij = —T^Uj^ — Tjk^i,k 
is the production term. Here, the Kolmogorov hypothesis of isotropy is invoked to model 
the dissipation rate tensor. 

The transport equation for the turbulent kinetic energy K = «/ 2 is obtained by 
contracting indices in Eq. (5) as 


DK 

Dt 


= V - £ + ltd 


( 6 ) 


where V = -rquq is the turbulence production. 
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As shown by Zeman and Sarkar, the dissipation rate of compressible homogeneous flow 
at high Reynolds numbers can be decomposed into a solenoidal or incompressible part, e s , 
and a compressible part, £ d as 

£ = £ s + £d (7) 

where e s = i/cj-ij- and £j = ) 2 given that V is the kinematic viscosity and lo\ is the 

fluctuating vorticity. Direct simulations of compressible, homogeneous shear flows show 
that e s is largely independent of compressibility and that the growth of turbulent kinetic 
energy decreases as the Mach number increases because of the augmented contribution of 
the compressible dissipation. Therefore, e d is modeled and e 3 is obtained by solving the 
incompressible form of the dissipation equation 


De s 

Dt 


^8 ~ 



( 8 ) 


where C sl and C e2 are closure coefficients that are model dependent. The models of the 
dilatational dissipation and pressure-dilatation correlation considered are 


Sarkar 


where M t = 
1.4). 


yJ^KhBf 


e c = 0.5 M?e s (9) 

p'd' = -0.1 5pVM t + 0.2 pe s M? (10) 

is the turbulent Mach number and 7 is the ratio of specific heats (= 


Zeman 


e c = 0.75(1 - txp(-((M t - 0.1)/0.6) 2 ))e, 


-1 ,P 


p'd' = (pry) X ( 


with 


T f 


—A) = 

Tf Dt v 7 p 


0.4— M t 

£> 


Pl = 2fK' f RT(— M t t M < ) 
1 1 +M? + M t 4> 


( 11 ) 

( 12 ) 

( 13 ) 

(14) 


Finally, we consider incompressible pressure-strain models, which are modeled as linear 
functions of the mean velocity gradients with coefficients that depend algebraically on the 
anisotropy tensor and the turbulent dissipation rate. Compressibility effects are incorporated 
through changes in density. Three models will be analyzed: the Launder, Reece, and Rodi 
(LRR) model (1975); the Fu, Launder, and Tselepidakis (FLT) model (1987); and the 
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Speziale, Sarkar, and Gatski (. SSG ) model (1991). These models are either linear (LRR), 
quadratic (SSG), or cubic (FLT) in the anisotropy tensor 

b - - ~ <15) 

and are assumed to be only functions of the anisotropy b and the symmetric and antisym- 
metric parts of the mean velocity gradient 


c 1 (foi | duA 

ij ~ 2 \dxj + dxj 

1 ( duj duA 

t; 2 \dxj dx, ) 

The high Reynolds number forms of these models are 


(16a) 

(166) 


Launder, Reece, and Rodi (LRR) 


n tJ (b,s,w) 


—Cieb tJ + C 2 K(Sij — \Skk&ij) + C3 A. (b l kSjk + bjkSik buSkAij, 


+ C4A (bikWjk + bjkWik) 

where 

Ci = 3, C 2 = 0.8, 

C 3 = 1.745, C 4 = 1.309 


(17a) 


(176) 


Fu, Launder, and Tselepidakis (FLT) 


n tJ (b,S,W) 


milVFe 


b z j + — ( b{kb, 




- bkibkAij ^ 


4 1 6 / 9 ~ 

+ —A ( Sij — — Skk ) + 7 A ( b t kSjk + bjkS t k — -bkiSkAij 

5 o 0 \ on 0 ' 

+ -A(6 lfc li^ + b jk W ik ) 
15 


+ -A' ^(bikbkiSji + bj k bk,S xl — 2b lk Subij — 36^5^6^) 

+ (bikbkiWji + bjkbkiWu)] 

— irK \siI(bikWjk +bjkWik) + 12 (bikbkiW trn b m j + 6 jt 6; c TT/ m 6m«)] 

L (18a) 
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where 


H 2 '■ III e^bijbjkbkl 

r = 0.7, F = 1 + 91 1 + 271 1 1 (18 6) 


Speziale, Sarkar, and Gatski (SSG) 


n tJ (b, S, W — [C\£ + C^V)bij + C 2 £ ( bikbkj — -bkibkiGj 


+ (c 3 - c; ] /b^b~) i<(Sij - 
+ C 4 K (bikSjk + bjkS’ik - -bkiSkibij'j 

+ C 5 K(bikWjk + bjkWik) 

where V = —TkiSki is the turbulence production and 


(19a) 


Cx = 3.4, C{ = 1.80, C 2 = 4.2, 

C 3 = 0 . 8 , c; = 1.30, C 4 = 1.25, ( 196 ) 

C 5 = 0.40 

As mentioned previously, the closure coefficients C £l ,C e2 , and C e for the turbulent dis- 
sipation rate equation (Eq. (8)) are model dependent. The values for these coefficients are 
given by 


Launder, Reece, and Rodi (LRR) 


C e i = 1.44, Ce 2 = 1-90 (20a) 

Fu, Launder, and Tselepidakis (FLT) 

C e \ = 1-45, C £ 2 = 1.90 (206) 

Speziale, Sarkar , and Gatski (SSG) 

Cel = 1.44, Ce2 = 1.83 (20c) 
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In the LRR and SSG models, the dissipation rate £ tJ is modeled through the usual isotropic 
assumption £ tJ = | whereas in the FLT model the tensor dissipation rate explicitly 
accounts for anisotropic effects: 

en = ley/FSij + 2(1 - V¥)eb xl ( 21 ) 

where F has been defined previously in Eq. (18b). 

Calculations of compressible, homogeneous shear flows show that the Reynolds stresses, 
the dissipation rate, and the mean temperature grow exponentially, so that the anisotropy 
tensor b,j, the shear parameter SK/s s and the turbulent Mach number M t achieve equilib- 
rium values that are independent of the initial conditions. Therefore, the system of equations 
for T x j , e, and T is nondimensionalized and is recast into an equivalent set of equations for 
e,/SK, b tJ , K, and M t as follows: 



m- 

(16) 


T' = m ‘<sX-< £ 

(17) 


dFs? = (^«(c..-<)-(^-( e e > 

(18) 

Dbij 

D t* 

~ 2SK f • e, 3e/’ J ^2 SI<^ ^ £ s £ s e 9 ^SK^ ,J 

(19) 


where t* — St is the dimensionless time, K* — K/Kq,Kq is the initial value of turbulent 
kinetic energy, and er = 7(7 — 1 )/2 . An additional equation, Eq. (12), is needed when 
Zeman’s model is used. 

This system of nonlinear ordinary differential equations associated with each Reynolds stress 
model is solved subject to the initial conditions 



0.45, M 2 t = 0.094, 


6 n = 0.124, b 22 = -0.106, b l2 = -0.187 


at time t* = 0, which is taken from the DNS of Blaisdell (run Shal92) at t* = 2. The reason 
for this is that the initial conditions of the DNS at f* = 0 are nonphysical. See Speziale 
et al. (1992). The system was integrated with a fourth-order Runge-Kutta scheme. The 
equilibrium states for each turbulence model are obtained numerically. A comparison of the 
predictions of the models with the DNS data will be made in the next section. 
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4. RESULTS AND DISCUSSION 


The predictions of the SSG, LRR, and FLT models with Sarkar’s compressibility correc- 
tions will be compared with the DNS results of Blaisdell (run Shal92). 

Figures 1 and 2 show the time evolution of the turbulent kinetic energy and the turbulent 
Mach number, predicted by the SSG model for two initial turbulent Mach numbers M t0 = 
0 and M t o — 0.307. These figures clearly show that a variable-density extension of the 
SSG model is not capable of capturing the decrease of the growth rate of turbulent kinetic 
energy as the turbulent Mach number increases. Explicit compressibility corrections are 
needed to predict the trends of the DNS. However, the model SSG overpredicts the growth 
rate of solenoidal and total dissipation (Figs. 3 and 4), which accounts for the overly large 
equilibrium value of the shear parameter (Fig. 5). Also, the differences in the predictions 
between turbulence models are shown (Figs. 6 - 10). The SSG and FLT models reproduce 
the results fairly well. However, the LRR model performs poorly because it was not well 
calibrated in incompressible, homogeneous shear flow. See Abid and Speziale (1993). 

In Figures 11 - 13, the model predictions for the time evolution of Reynolds stress 
anisotropies are shown. All of the models substantially underpredict the Reynolds stress 
anisotropies in comparison with the results of the DNS. Also, the relaxation of b t j to their 
equilibrium states has not been captured by all the models. As pointed out by Speziale, 
Gatski, and Sarkar (1992), the long time behavior of turbulence models is tied to their 
ability to predict the equilibrium values. 

In fact, all the models predict exponential long time behavior, i.e., K, e, T, and 
proportional to exp(X 00 t m ), where A^ is the equilibrium growth rate given by 

a » = -2( m » - C-i-jr)* ( 2 °) 

The equilibrium states obtained from the various turbulence models are compared with 
the DNS data of Rogers et al. (1986) for incompressible homogeneous shear flow (table 
I) and the DNS data of Blaisdell (1990) for compressible homogeneous shear flow with 
Sarkar’s model (table II). Several observations in regard to these results are noteworthy: the 
Reynolds stress anisotropies (6^)^ and the shear parameter [SK/e s ) oo are underpredicted 
by all models, particularly (6n)oo; the shear stress (612)00 is erroneously predicted to be 
insensitive to the compressibility effects; and the SSG and LRR models predict the observed 
trend that (6n)oo increases as a function of the turbulent Mach number M t (however the 
FLT model erroneously predicts the opposite trend of the DNS data). Similar conclusions 
are drawn when Zeman’s model is used. See table III. 

All of the turbulence models are clearly incapable of properly accounting for the effects 
of compressibility on turbulent shear flow. The failure of the Reynolds stresses to predict 
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the increase in the anisotropy of normal stresses and the reduction in the shear stress is 
partly due to the use of the incompressible pressure-strain correlation models. As shown by 
Blaisdell and Sarkar (1993), the contribution of the compressibility to the pressure-strain 
correlation is large and must be taken into account. 

Another deficiency of the turbulence models considered in this study lies in the use of 
the incompressible form of the dissipation equation to obtain the solenoidal part of the 
dissipation rate. The results presented in tables I and II show clearly that ("P /£ s )oo is 
erroneously predicted to be sensitive to the compressibility effects. When an equilibrium 
state is achieved, Eq. (18) gives 



C e 2-1 

C.1-1 


Cel-1 


)( 


£d ~ Kd, 


( 21 ) 


where the first term in the right side represents the equilibrium value of V /e s in incompress- 
ible, homogeneous shear flow. Equation (21) clearly shows that the erroneous prediction 
of {V/£ s )oo is directly tied to the solenoidal dissipation equation. This conclusion can be 
reached in other ways. 

The value of (P/e s )oo can be determined from Eq. (17) as 


(-)oo = (l + ( £i - J %l + -7AO (22) 


The combination of Eqs. (21) and (22) leads to 





£2 


( c <i - 1 + n2w>”) 


(23) 


The above equation, which is independent of the dilatational dissipation and pressure- 
dilatation models, shows that (P/e s )oo cannot be insensitive to the compressibility effects 
as seen in the direct simulations, otherwise Eq. (23) will predict nonphysical values for M t 
(M too = 1.7 when ( V/£ s )oo — 1-84). Note that the system of Eqs. (16)-( 19) will predict 
nonphysical results for M too if the explicit compressibility equations are not included. 


4. CONCLUDING REMARKS 

Three Reynolds stress models have been evaluated for the problem of compressible, homo- 
geneous shear flow. The dilatational dissipation and pressure-dilatation correlation models, 
formulated by Zeman and Sarkar, have been inserted in the Reynolds stress closure. Com- 
parisons between the predictions of the various models and the direct numerical simulation 
of Blaisdell have been made. 

All three Reynolds stress models fail to capture the compressibility effects seen in the 
direct simulations. In particular, the increase in anisotropy of the normal stresses and the 
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reduction in the shear stress are not reproduced by any of the models. This result is partly 
due to the use of variable density extensions of incompressible pressure-strain models. Thus, 
dilatational effects on the pressure-strain correlation must be identified and accounted for in 
compressible turbulence modeling. 

An analysis of the predictions of the equilibrium states has shown that use of the in- 
compressible form of the dissipation to determine the solenoidal dissipation rate is another 
source of inaccuracy in the predictions. In particular, the erroneous prediction of (P/e 5 )oo 
is directly tied to the solenoidal dissipation equation. 

The present study indicates the need for further direct numerical simulations of com- 
pressible, homogeneous shear flows. Such simulations could be used to distinguish between 
compressibility and low Reynolds number effects and to provide more information on the 
equilibrium states. 


10 


REFERENCES 


Abid, R. and Speziale, C. G., 1993, “Predicting equilibrium states with Reynolds stress 
closures in channel flow and homogeneous shear flow,” Phys. Fluids A 5, pp. 1151- 
1158. 

Blaisdell, G. A., 1990, “Numerical simulations of compressible homogeneous turbulence,” 
Ph.D. Thesis , Mechanical Engineering Dept., Stanford University, Stanford, California. 

Blaisdell, G. A. and Sarkar, S., 1993, “Investigation of the pressure strain correlation in 
compressible homogeneous turbulent shear flow,” in FED 151, Transitional and Tur- 
bulent Compressible Flows. 

Fu, S., Launder, B. E. and Tselepedikas, D. P., 1987, “Accommodating the effect of high 
strain rates in modeling the pressure-strain correlation,” UMIST Mechanical Engineer- 
ing Department Report TFD/87/5. 

Launder, B. E., Reece, G. and Rodi, W., 1975, “Progress in the development of a Reynolds 
stress turbulence closure,” J. Fluid Mech. 68, p. 637. 

Rogers, M. M., Moin, P. and Reynolds, W. C., 1986, “The structure and modeling of 
the hydrodynamic and passive scalar fields in homogeneous turbulent shear flows,” 
Technical Report TF-25, Stanford university, Stanford, California. 

Sarkar, S., Erlebacher, G. and Hussaini, M. Y., 1991, “Direct simulations of compressible 
turbulence in a shear flow,” Theoret. Comput. Fluid Dynamics 2, pp. 291-305. 

Sarkar, S., 1992, “The pressure-dilatation correlation in compressible flows,” Phys. Fluids 
A 4, pp. 1251-1258. 

Speziale, C. G., Sarkar, S. and Gatski, T. B., 1991, “Modeling the pressure-strain correlation 
of turbulence: An invariant dynamical systems approach,” J. Fluid Mech. 227, p. 245. 


Speziale, G. S., Gatski, T. B. and Sarkar, S., 1992, “On testing models for the pressure-strain 
correlation of turbulence using direct simulations,” Phys. Fluids A 4, pp. 2887-2899. 

Zeman, 0., 1990, “Dilatation dissipation: The concept and application in modeling com- 
pressible mixing layers, ” Phys. Fluids A 2, p. 178. 

Zeman, 0. and Blaisdell, G.A., 1991, “New physics and models for compressible turbulent 
flows,” in Advances in Turbulence 3, A. V. Johanson and P. H. Alfredson eds., Springer 
Verlag, Berlin, pp. 445-454. 


11 



Equilibrium 

Values 

LRR 

Model 

SSG 

Model 

FLT 

Model 

DNS 

Data 

bn 

0.155 

0.219 

0.208 

0.215 

b 12 

-0.187 

-0.164 

-0.146 

-0.158 

b 22 

-0.121 

-0.146 

-0.144 

-0.153 

b 33 

-0.004 

-0.073 

-0.064 

-0.062 

SK/e. 

5.34 

5.77 

6.84 

5.70 

P/e s 

2.0 

1.88 

2.0 

1.80 


Table 1. Predicted Equilibrium Values for Incompressible Homogeneous Shear Flow. 
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Equilibrium 

Values 

LRR 

Model 

SSG 

Model 

FLT 

Model 

DNS 

Data 

611 

0.166 

0.230 

0.189 

0.424 

612 



-0.148 

-0.118 

622 



-0.138 


&33 

-0.036 

-0.082 

-0.051 

-0.188 

SK/e, 

3.77 

4.11 

4.77 

7.82 

M t 

0.65 

0.60 

0.65 

0.51 

P/e s 

1.41 

1.36 

1.41 

1.84 


Table 2. Predicted Equilibrium Values for Compressible Homogeneous Shear Flow with 
Sarkar’s Model. 
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LRR 

Model 

SSG 

Model 

FLT 

Model 

DNS 

Data 

0.167 

0.231 

0.187 

0.424 

-0.191 

-0.167 

-0.148 

-0.118 

-0.131 

-0.148 

-0.137 

-0.236 

-0.036 

-0.083 

-0.050 

-0.188 

3.59 

3.95 

4.61 

7.82 

0.48 

0.45 

0.48 

0.51 

1.41 

1.36 

1.41 

1.84 


Compressible Homogeneous Shear Flow with 
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Figure 1. Comparison of the predictions of the SSG model for the time evolution of the 
turbulent kinetic energy with the DNS results of Blaisdell. 
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Figure 3. Comparison of the predictions of the SSG model for the time evolution of the 
solenoidal dissipation-rate with the DNS results of Blaisdell. 
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Figure 4. Comparison of the predictions of the SSG model for the time evolution of the total 
dissipation-rate with the DNS results of Blaisdell. 
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Figure 8. Comparison of the predictions of the LRR, SSG and FLT models for the time 
evolution of the solenoidal dissipation-rate with the DNS results of Blaisdell. 
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Figure 9. Comparison of the predictions of the LRR, SSG and FLT models for the time 
evolution of the total dissipation-rate with the DNS results of Blaisdell. 
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Figure 12. Comparison of the predictions of the LRR, SSG and FLT models for the time 
evolution of the component b n of the Reynolds stress anisotropy tensor with the DNS results 

of Blaisdell. 
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